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Abstract 

The dynamics of amphiphilic films deposited on a solid surface is analyzed 
for the case when shear oscillations of the solid surface are excited. The two 
cases of surface- and bulk shear waves are studied with the film exposed to 
gas or to a liquid. By solving the corresponding dispersion equation and the 
wave equation while maintaining the energy balance we are able to connect 
the surface density and the shear viscosity of a fluid amphiphilic overlayer 
with experimentally accessible damping coefficient, phase velocity, dissipation 
factor and resonant frequency shifts of shear waves. 
PACS number(s): 68.15.+p, 68.60.Bs, 43.35.Bf 
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I. INTRODUCTION. 



Measurements of the properties of thin adsorbed films have since long been an important 
subfield of interface studies [] 10 1 . Scientific and technological applications of thin films like 



Langmuir-Blodgett (LB) films [ ! L — 2C ] , self-assembled monolayers (SAM) |21[] and protein 
monolayers and multilayers |2"2"|J2"3"1 have stimulated detailed studies of the physical prop- 
erties of such structures. With these model systems it is possible to construct multilayer 
amphiphilic films with controlled monolayer thickness hj ~ 25A and to create different cou- 
pling situations between the adsorbed film and the chemically modified substrate |T2}-|To|. 
Specifically measurements of surface density and viscosity of these films are important both 
scientifically and, e.g., with respect to their possible application in acoustical biosensors in 
which they can serve as sensitive elements jl3|,|l6|Jl8| , |l9f . 

The dynamic (shear) viscosity of amphiphilic films is a very important rheological char- 
acteristic of a protein film or a fluid membrane, which is strongly dependent on temperature 
and phase state of the film, and on environmental conditions |T^j2^j25[| . Usually LB or 
SAM film acoustic sensors operate in vacuum or in a gaseous enviroment. In contrast, a 
lipid-water system or a protein layer adsorbed onto a solid surface in a fluid experiences 
"wet" conditions or a bulk aqueous medium [23,26] at room temperature, which can change 
the viscosity of the film. In addition, the biomolecular layer adsorption process is extremely 
sensitive to the nature of the solid substrate in the sense that it can modify the overlayer 
structure HJUJiMII . In 

view of this, a rigorous description of the dynamics of the existing 
interfaces of combined solid substrate-adsorbed layer-bulk liquid systems, as well as of solid 
surface-amphiphilic film-gas systems seems to be greatly needed. 

In the present paper we investigate the dynamics of thin amphiphilic layer(s) attached 
to a solid substrate oscillating in a shear mode and in contact with (i) a gas or (ii) a bulk 
newtonian liquid. Figs. (1-3) show schematically the geometry of such sandwich systems. 

Amphiphilic bilayer films can imitate the behavior of a bilayer lipid membrane. The 
constituent bilayer molecules are composed of hydrophilic groups attached to hydropho- 
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bic chains of different length. The phase diagrams of such amphiphilic molecules demon- 
strate a rich variety of properties and behavior depending on temperature and water content 
[13,20,24,27,28]. Above the liquid crystalline-to-gel transition the hydrocarbon chains are 
approximately liquid and a bilayer membrane behaves like a two-dimensional fluid in the 
lateral plane due to the vanishing shear modulus of elasticity [27,28]. Such a peculiar fluid 
is isotropic in its plane but it is anisotropic in the normal direction due to sublayered head- 
and-tail structure [24]. Viscous layer response on shear deformation can be characterized 
by different shear viscosity coefficients, namely, the surface viscosity, n s [13], and a bulk 
shear viscosity, t)m (a component of tensor rjjkim, index 'M' corresponds to the classification 
according to Ref. 28). 

When the system is immersed in water, anisotropy can arise from water being trapped 
in the layer. This can result in a value of the overlayer viscosity n M , which is distinctly 
different compared to the case of a gaseous (air) enviroment. This is especially important 
for protein adsorbed layers with nonuniform interfacial domain structure and nonuniform 
water distribution in the normal direction. Such nonuniformity effects have been analyzed 
in a recently published hydrodynamic model of "porous" polymer surface film in liquids (for 
a review see Ref. 26 and references therein). 

In the present paper we have considered both the surface viscosity n s and the bulk shear 
viscosity i]m of adsorbed amphiphilic films. These can be measured in different acoustic 
experiment, when acoustic shear waves propagate along the plane of the layer (horizontally 
polarized surface shear waves) and in the direction normal to it (bulk acoustic shear waves). 
In order to be able to describe both air- and liquid enviromental experiments, corresponding 
shear viscosity coefficients, n M 'in air' and n M 'in liquid! have been noted. 

The paper is organized as follows: in the framework of continuum mechanics we derive in 
Section II the dispersion equation for elastic shear waves propagating in the system. We also 
treat the interaction between horizontally polarized elastic surface shear waves (SH SSW) 
and a fluid layer adsorbed on the surface of a semi-infinite solid substrate. In Section III 
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we analyze the response of bulk acoustic shear waves (BAW) propagating in a finite quartz 
plate with a viscous overlayer. Our results allow us to connect the surface density and 
shear viscosity of amphiphilic overlayers with experimentally accessible damping coefficient, 
phase velocity, dissipation factor and resonant frequency shift of shear waves. These can be 
measured by modern piezo-acoustical devices of different types [^,^,^]. 



II. DYNAMICS OF A VISCOUS BILAYER FILM ON A SOLID SUBSTRATE 
OSCILLATING IN A GASEOUS MEDIUM OR IN VACUUM. 

Within continuum mechanics [30,31] the dynamics of a liquid film can be described by 
the equation 

dv ' 

p—± = p (dvj/dt + (vV) vj) = d k a jk . 

This corresponds to the Navier-Stokes equation for the motion of a viscous liquid. Recently it 
was shown that hydrodynamical properties of fluid amphiphilic (lipid) films can be analyzed 
within the continuum mechanics scheme with viscous tensions defined as follows [28]: 

Aim = \ (div m + d m vi) , di = d/ dxi . 
Here rj^im is a viscosity matrix which were introduced explicitly in Ref. [[52 



For a fluid lipid layer in the incompressible liquid approximation and in the absence of 
pressure gradient, these expressions can be simplified p8 |. In particular, for the case when 
a strain is applying along the x direction, one finds 

Vxy = VM (d y V x ) . 

Here the viscosity coefficient t]m is a component of the viscosity matrix %Mmi the y-axis is 
perpendicular to the layer plane. 

In our model we use a notation t]m (liquid) and T]m (gas) for describing fluid films 
bounded (in different experiments) by a liquid or a gas. 



Below we analyze the influence of viscosity and surface density of a thin bilayer on the 
phase velocity and the damping of elastic shear waves. 



Surface Shear Wave Propagation: The Semi-infinite Quartz Crystal Substrate 

We treat here horizontally polarized surface elastic waves on a solid substrate interacting 
with a thin fluid double layer, which on the other side has an interface to vacuum or to a 
gaseous phase (Fig.l). Surface waves on plane interfaces have an amplitude which decays 
exponentially with normal distance from the solid surface on which they propagate [2]. Phase 
velocities of surface shear waves (SSW) V are lower than those of bulk shear waves Vo in 
the semi-infinite elastic half-space (i.e. the substrate). 

Let us now consider how the shear vibration of the substrate generates a viscous wave 
in the adjacent fluid layer 0,[|,|17]]. In the framework of fluid mechanics, the motion of an 



adsorbed viscous film is described by the linearized Navier-Stokes equation: 

dv x (y,z,t)) jdt = viAv x , 

^i = Vi/Pi, W 
Vi = Vm (gas). 

Equation (1) is valid in the regime of small Reynolds numbers. 

R = ^h <<1; 

V\ 

for small oscillation frequency u and -amplitude uq and for a thin overlayer thickness h\. Here 
v x denotes the x-component of the velocity of the fluid film and the ?/-axis is perpendicular 
to the z-direction of wave propagation. 

The boundary conditions at the fluid-solid interface (y = 0) correspond to the assumption 
of no slip conditions [6,26,28,30] 

y = o : Vx = 

vv=m%, ( 2 ) 

_ (~t du x 
Vyx — ^44 dy ■ 



Here u x is a component of the substrate boundary displacement vector, 

u x = u exp(K t y + iqz)exp(iujt), 

and Kt = (q 2 — pou> 2 /C^) 1 ^ 2 is the inverse penetration depth of surface shear waves into the 
solid substrate; u, q are wave frequency and wave number, while po and C44 are the substrate 
density and shear modulus, respectively. The substrate boundary motion is described by 
the equation 

,d 2 u x d 2 u x . 
p u x = c 44(^- + -Q^-h 

which corresponds to elastic vibrations of the solid surface. 

Boundary conditions at the moving interface between fluid 1 and fluid 2 (y = h±, Fig. 1) 
follow from the condition that the friction forces must be the same [[3(]] : 



U C 1 ) (2) 

y = h u n k a y ik ' = n k a\ k ' , 

W - A., J_ f dv^l j_ dvklll 

ik 



(3) 



<7jfc - -P2<>ik + V2{^- + "fe-J, 

w w = „« = „w = 0, 

At the free surface, finally, the boundary condition is: 

fa® 

V = ha, ifc-^- = 0, (4) 

Equations (1-4) together with the equation for the elastic substrate motion lead to the 

following dispersion equation for horizontally polarized surface shear waves in our system: 

iurqib (| + etanh(A^ 2 ))exp(2^i) - (| - etanh(Ah&)) 
Kt ~ C u (| + etanh(A^ 2 ))exp(2^ 1 ) + (|-etanh(A^ 2 ))' ^ 

where £j = (q 2 + icu/uj) 1 ^ 2 , j = 1, 2, e = 772/771 , and A/t = h 2 — h\. 

Here we introduce the viscous penetration depth 5 = (2z//c ( j) 1 ^ 2 corresponding to the 
distance over which the transverse wave amplitude falls off by factor of e. For two sufficiently 
thin viscous overlayers one can assume that hj/5j <C 1, and we will always consider the long 
wavelength limit, where q 2 <C 5~ 2 . Within these limits, we can from Eqn. (5) find the SSW 
damping coefficient as a function of the viscosities and the surface densities of the overlayers. 
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The result for the damping coefficient T of the SSW — simply the imaginary part of the 
wave vector q — is 

r = r^{i + *(i + ^(i + *YU. ( 6) 

{ Vis \ V2spls V VuJJ) 

Here T]j S = r] M hj is the surface viscosity component of each layer, pj s = pjhj is a correspond- 
ing surface density and T mon denotes the damping coefficient for the case of a monolayer 
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q uj 3 p ls n ls 

-5 



One can estimate the value of T to be ~ 10 for a phospholipid bilayer. Experimental data 
for egg lecithin at T = 25° C are taken from Ref. ||24|| , which provides the following values: 



bilayer thickness hb = 46A, bilayer density p ~ 1 g/cm 3 , bilayer "microviscosity" (measured 
by a probe technique) n « 1.2 dyn-s/cm 2 . The resonator frequency is c^o ~ 10 9 Hz and 
C44 ~ 2.6 dyn/cm 2 is the corresponding quartz shear modulus. 

Using the dispersion equation (5) one can also obtain the change in the SSW phase 
velocity, AV/V , caused by the presence of the adsorbed bilayer. One finds that 

^„ + *-,» (7) 

V 2 Vopo p 3l 



Vb = \JC 44 /p = const. 

As a concequence of the no slip- and thin layer assumptions, the SSW velocity shift 
is sensitive only to the overlayer surface densities p Sj and not to their viscosities. This 
corresponds to the Love type of wave propagation. Equation (7) allows one to determine 
the surface density of the upper or lower half of the bilayer if the density of the other layer is 
known or can be determined by an independent experiment. For the lipid bilayer (with the 
same parameters as above) the SSW phase velocity shift is small, AV/Vq ~ 10~ 8 . However, 
this shift may be detectable (of order ~ 10~ 6 ) at high frequences, u ~ 27rl0 9 Hz, since the 
SSW velocity change is a quadratic function of frequency. Such high frequencies can be 
generated by modern acoustic SSW devices |29|. At low frequences, /^10 8 Hz, it is more 
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reasonable to use another type of shear waves for surface density analysis excited by quartz 
crystals oscillating in the thickness shear mode (see Section III.B below). 



III. SOLID SUBSTRATE WITH VISCOUS OVERLAYER OSCILLATING IN A 

BULK VISCOUS MEDIUM. 

Resonant Frequency and Dissipation Factor for Bulk Shear Waves: The Finite 

Quartz-Crystal Resonator 

In addition to "genuine" surface waves, another type of shear waves can propagate in the 
plane of the overlayer-substrate interface. These are bulk acoustic waves (BAW). In contrast 
to the SSW case considered above, the acoustical response of unloaded BAW resonators 
depends on the quartz slab thickness h as well as on its density p and shear modulus C 66 

fo = (C 66 /po) 1/2 ^o 

It has been shown []^,|6l p3l,p6l|3~3l that the resonant frequency of such quartz plates de- 



creases when its surface is coated with an overlayer constituting a mass load. For a thin 
surface film, uniformly covering the entire vibrating area of the slab, Sauerbrey deduced a 
linear relation between the frequency change and the added mass per unit area in a vac- 
uum (gaseous) environment [Q . Kanazawa and Gordon || have later described the resonant 
response of quartz oscillators in bulk viscous liquids. Recent experimental and theoreti- 
cal investigations of the BAW in composite media (see Refs. 23, 26, 33 for a review) have 
demonstrated possible applications to the study of nonuniform polymer and protein thin 
films and have included both resonant frequency and viscous dissipation in the interface 
region. Experimentally, it is convenient to obtain this viscous dissipation as the width of 
the resonant frequency [26] or as the dissipation factor D [23] 

D = Edissipated/ZnE stored = I/ 71 "/ 7 " 

where r is the time constant for the decay of the vibration amplitude. 



In this Section we calculate the shifts Af of the resonance frequency and AD of the 
dissipation factor, respectively, for the case of a viscous amphiphilic (lipid) film attached 
without slip to the surface of the quartz slab oscillating in the thickness shear mode in the 
presence of a bulk viscous liquid. 

The wave equation for shear waves propagating in the vertical direction (Fig. 2) is 

d 2 u x (y,t) ipmj 



dy 2 

Its general solution has the form 



-u x (y,t), 771 = T] M (liquid) . 



(8) 



u , 



e iuJt (U^v + U 2 e^- 

ei = (i+i)M 



(9) 



Using the same boundary conditions (2-4) as above, we find for the x-component of the 
liquid velocity: 

e Ziy _|_ AQ-Ziiv-^hi) 



v x = v - 



1 + Ae 2h ^ 



(10) 



62 + eSi tanh 


+ Ah/62 


82 — eS\ tanh 


+ Ah/62 



A 



The shift of resonant frequency Af and dissipation factor AD, due to the overlayer on 
the substrate, can be calculated from the balance between dissipated and stored energy in 



the system [33,26,23]. Using this energy balance, we find from Eq. (11): 



Af « -Im 



Ae 2h ^ + 1 



/2irp h 



(11) 



/ Ae 2hl ^-l\ 
AD « -Re Vl ^——-— /nfp h 



Ae 2h ^ + l r 

In the limiting case h\/8\ <C 1, Ah/ 82 > 1 we obtain for the resonant frequency shift 



(12) 



, / fV2p2 
\TTp C 66 



1/2 



1 + /lipi 



' knV, -(*r/^' 



JI2P2 



m 



(13) 



This expression contains contributions from both the lower and the upper layer. For suf- 
ficiently viscous solutions, r] 2 ~ rjipi/ p 2 , the second and third terms may be of the same 
order. 

In contrast, the result for AD is practically insensitive to the ultrathin overlayer param- 
eters. In the linear approximation, corresponding to a small value of hi/ Si, it depends only 



on the upper (bulk) liquid viscosity and its density: 

1/2 

ii / ; iii i /;/.,/ i.-i \ i 

(14) 



\^JPoC m ) { rji \rjipi , 



In the limit of r^ 2 — > we obtain the Sauerbrey formula: 

A/=-^. (15) 

V ^0^66 

In the opposite case and for hi = 0, the results of Kanazawa and Gordon follow: 

v-A&T- (16) 

Formulae (14) and (15) form a set of equations. From these expressions it is possible to 
determine the i] M component of the overlayer shear viscosity and the thickness (mass) of the 
adsorbed layer and compare them with the results of surface viscosity measurements. Also, 
it is possible to determine the viscosity and density of the adjacent bulk liquid if parameters 
of the overlayer are known. 

It is essential to note that our model is valid for an arbitrary ratio e between the viscosities 
i]2 and r]i, i.e., even for e > 1. The latter possibility may be relevant for a thin amphiphilic 
substrate-adjacent film under a thick protein layer (For example, a sandwich structure such 
as the one schematically shown in Fig.4. can be realised by an adsorbed protein layer 
attached to the substrate via an amphiphilic self-assembling monolayer [22]). 



IV. DISCUSSION. 



In addition to its importance in technological applications, the study of the hydrody- 
namic modes of fluid amphiphilic films and adsorbed proteins are generally important for 
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understanding the dynamics of the swelling of lipid-water systems [28,34,35] and the rheol- 
ogy of biological membranes. It is known, for instance, that living cell membranes, e.g. in 
red blood cells, experience shear stress in hydrodynamic flow through capillaries [^,0. As 
the lipid bilayer matrix of the membrane is a two-dimensional incompressible liquid adjacent 
to an elastic protein (spectrin-actin) network, the hydrodynamics of this layered structure 
is governed by the coupling of the fluid membrane to shear flows in an external bulk liquid 



The surface viscosity of a fluid amphiphilic film is often considered to be a two- 
dimensional analogue of the bulk viscosity and defined as the coefficient of proportional- 
ity between the tangential force per unit length and the gradient in the flow velocity of 
the liquid. These averaged surface characteristics can be measured, e.g. by an interface 



shear rheometer at air-liquid interfaces [|T|,[25[ or by using the oscillating barrier method in 
a LB trough [25]. In contrast to r] s , the "microscopic" shear viscosity may be defined as 
a contribution of fluid disordered chains. The "microviscosity" of the membrane core can 
be determined experimentally, e.g. by a probe technique, and describes the local viscos- 
ity near the probe |23j. This "microviscosity" is often identified with the true membrane 
viscosity. However, this experimental technique could give rather different results with a 
strong dependence on the probe material and on the interaction between the probe and its 
local surroundings f24|. Both of these viscosities may be changed after the deposition of an 
amphiphilic layer on the solid surface. 

The results of our theory can be applied for direct acoustic measurements of surface 
densities and/or shear viscosities i] s and t)m for a thin amphiphilic film after the transferring 
onto the solid substrate and in both the liquid and gaseous experimental conditions. Another 
help is a possible strong coupling between the end groups of lipid layer, LB film or SAM 
and the surface of the substrate [12-16], thus for the interfacial solid-fluid region, no slip 
assumption seems to be valid. 

In our work we found the dispersion equation for surface shear acoustic waves with 
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horizontal polarization and solved it as well as the wave equation for shear bulk acoustic 
waves. The solution of these equations together with energy balance, allows us to calculate 
the analytical expressions for damping coefficients, phase velocities, dissipation factor and 
resonant frequency shifts of both types of shear waves as functions of lipid film surface 
density and two different components of its in-plane viscosity. These acoustic waves can be 
excited by means of two different types of piezoelectric oscillators [23,29]. 

We can suggest to measure the acoustical responce of the above-mentioned sandwich 
system in the region of lipid phase transition when the viscosity of the overlayer subject to 
dramatic changes. For instance, during the liquid crystalline to gel transformation of a lipid 
membrane, the viscosity of the bilayer changes by more than 10 times while the membrane 
density remains practically constant [17,24]. It is an interesting experimental fact [34] that 
the two halves of a bilayer lipid membrane are so weakly coupled, that they can udergo 
the thermotropic phase transition independently. In the accordance with our results, it 
must provide the changes in SSW damping coefficient, while the corresponding SSW phase 
velocity shift due to the presence of bilayer will be constant. 

On a more speculative note, our results may also be valid within the segregated (sublayer) 
structure model of adsorbed protein layers (Fig. 3). In accordance with this model (see, 
e.g. [21,22]), the structure of the protein overlayer resembles a surfactant film composed 
of segregated head-and-tail regions. Recent neutron scattering experiments |22|] have given 
evidence for such an anisotropic sublayer structure of protein monolayers adsorbed onto solid 
surfaces through the self-assembling amphiphilic sublayer. In this case, the relation between 
the upper and lower layer viscosities may vary in strong dependence of the phase state of 
protein and amphiphilic monolayer, respectively. 
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FIGURES 




FIG. 1. Illustration of surface shear waves (SSW) with horizontal polarization propagating 
in a system of an elastic half-space with an adsorbed thin amphiphilic fluid bilayer on top. The 
system is in a gaseous enviroment. 



17 



A 



-h 



bulk liquid 



ox) o o o o o o o o°o9 
n ao n aQ.o.ao n o n (iQ.Q() : 

o%Wo o o o o o o ao flt 
o 








FIG. 2. shear bulk waves propagating in an AT-cut quartz plate covered with a thin am- 
phiphilic layer. The system is immersed in a bulk liquid. 
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FIG. 3. Schematic depiction of two distinct regions of protein layer absorbed from bulk 
solution onto a solid substrate. The dotted line correspondes to the surface-adjacent amphiphilic 
layer. 
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